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Quantum states are successfully reconstructed using the 
maximum likelihood estimation on the subspace where the 
measured projectors reproduce the identity operator. Recon- 
struction corresponds to normalization of incompatible ob- 
servations. The proposed approach handles the noisy data 
corresponding to realistic incomplete observation with finite 
resolution. 



Quantum theory handles observable events on the 
most fundamental level currently available predicting the 
statistics of quantum phenomena. This randomness is 
hidden in the quantum state. Although the history of 
state reconstruction may be traced back to the early days 
of quantum mechanics to the Pauli problem Q , till quan- 
tum optics opened the new era. Theoretical prediction of 
Vogel and Risken 1^ was closely followed by the experi- 
mental realization of the suggested algorithm by Smithey 
et. al. 1^. Since that time many improvements and new 
techniques have been proposed, an overview can be found 
in Ref. Q. Reconstruction of quantum states is consid- 
ered as standard technique used in various branches of 
contemporary physics [p|. Unfortunately, there are sev- 
eral flaws in the approach. The available measurement is 
always limited as far as the amount and accuracy of data 
is concerned. Standard methods are designed for analy- 
sis of sharp, complete and noiseless observations. When 
applied on the realistic data, serious problems with pos- 
itivity of reconstructed density matrix appear. Though 
these techniques may give a rough picture of the state, 
they cannot provide a full quantum description. 

This may be accomplished using the maximum like- 
lihood (MaxLik) estimation The question of deter- 
ministic schemes: "What quantum state is determined 
by that measurement?" is replaced by the formulation 
consistent with quantum theory: "What quantum states 
seem to be most likely for that measurement?" Quan- 
tum theory predicts the statistics provided that quantum 
state is known. The fundamental result of this Letter is 
the statistical inversion of this quantum postulate pre- 
dicting the quantum state provided that results of the 
measurement are known. General theory is formulated 
for the case of nonorthogonal measurement of incompat- 
ible observables. Physically it corresponds to synthesis of 
various measurements done under different experimental 
conditions, always performed on the identically prepared 
system. It might be subsequent recording of an unknown 



spin of the neutron using different settings of the Stern 
Gerlach apparatus, or the recording of the quadrature 
operator for light in the rotated frame in quantum to- 
mography. The formulation provided here will focus the 
general aspects illustrating the novel approach on the ex- 
ample of quantum tomography. 

In the following the states \yi) will denote general 
nonorthogonal states enumerated for concreteness but 
without loss of generality by a discrete quantum num- 
ber. Assume that the quantum measurement has been 
done n times yielding the relative frequencies of the 
event represented by the projector \yi) as fi > 0. The 
states are assumed to be nonorthogonal and linearly in- 
dependent. Hence the correlation matrix is invertible 
and Hermitian dj = C*, = {yi\yj)- (i) The measure- 
ment is sharp provided that it may be described by a 
projector into the pure state \yi){yi\. On the contrary, 
projectors corresponding to unsharp measurement are 
given by the probability operator measure [Q| represent- 
ing a superposition in bins Di of indistinguishable states 
Ilj — J2j<£D \yj){yj\- i^"^) Measurement is complete if all 
the projectors corresponding to the counted data yield 
the decomposition of unity operator, X^i a;; \yi){yi\ = 1- 
The states are assumed to be normalized with respect 
to this full completeness relation but not to the scalar 
product as usually. Provided that the counted data do 
not exhaust all the values, the measurement is incom- 
plete. In the notation used here, the decompositions are 
assumed to be incomplete unless defined otherwise, (in) 
The data are noise free provided that the counted fre- 
quencies coincide exactly with the prediction of quantum 
theory 



f^ pu = {y^\p\yi) 



(1) 



for some density matrix p and data are noisy otherwise. 
Sharp, complete and noise free measurement is assumed 
implicitly in the standard description of quantum to- 
mography based on the direct inversion of ([l]). When 
this ideal conditions are not met, the algorithm does not 
guarantee the positive definiteness of the reconstructed 
"density matrix" . 

Realistic data can never provide complete information 
about quantum system with infinite degrees of freedom. 
This may be demonstrated on the simplified example of 
reconstruction of diagonal elements of density matrix via 
photon counting with ideal detector. Suppose n times 
repeated counting, always with the zero registered pho- 
toelectrons. The "standard" prediction of quantum state 
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reads p = |0)(0|, where |0) denotes the vacuum state. 
Nevertheless this interpretation may be always spoiled 
by the classical noise represented by a projector into 
the strong coherent state = |-^)co/i(-^| appearing 
with the negligible probability e. Obviously, the state 
/5e = (1 — e)/5 + eAAc cannot be distinguished from the 
standard one for sufficiently small e < 1/n. This may 
appear as crucial for some observations. For example, 
the average numbers of particles differ significantly for 
both the states. The state of the system is well defined 
as the "standard" density matrix on the scanned part 
of the Hilbert space spanned by the actually measured 
values. The behavior of the wave function in the comple- 
mentary world is t unknown. If some future detection of 
a quantum variable will depend significantly just on this 
"unobserved" part, the prediction should be uncertain. 
This seems to be evident in the case of orthogonal mea- 
surements but as will be seen, rather nontrivial in case 
of nonorthogonal measurements. 

Quantum state attached to the data will be searched 
as the density matrix p which maximizes the likelihood 
functional Q 



C{p) ^X{{y^p\yr^^ 



(2) 



Assume the diagonal representation of a density matrix 



^^rk\(l)k){<f>k\- 
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The existence of parameters > 0, = I and an 

orthogonal basis is guaranteed by quantum theory. 
Normalized extremum states satisfy the relation 
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A being a Lagrange multiplier. This reads the system of 
coupled equations 



(4) 
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Pii = ^'rk\{(f)k\yi 



In the derivation, the condition of normalization Tr/5 = 1 
has been used. Relation (^) provides the statistical inver- 
sion of quantum postulate (l|) as the nonlinear equation 
for density matrix 



R{p)p ^ p. 



(5) 



Reconstruction is done in the subspace where operator 
R represents the identity operator. The equations for 
matrix elements pij read 



fiPi] = Pii^C^^^pkj, 

k 



(6) 



C ^ being the inversion matrix to Cy ={Vi\yj)- Notice 

i) 



that the diagonal elements are instead of 
relation 



y) fulfilling the 
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Though the relation (|^) is linear with respect to the den- 
sity matrix, the inversion represented by (^) or by (^) is 
not. The reasons are fundamental: Elements of density 
matrix are not independent, but characterize a quantum 
state. These quantum correlations are neglected when 
the inversion is done regardless on the positive definite- 
ness. To find solutions of nonlinear operator equation 
is a peculiar problem. It may be approached itera- 
tively provided that necessary conditions for convergence 
are fulfilled. This questions will be addressed separately 
elsewhere. Formulation simplifies considerably provided 
that the projectors in operator R commute. The density 
matrix is diagonal in this common basis. Reconstruction 
of diagonal elements of density matrix using the homo- 
dyne detection with random phase represents an explicit 
example ||^ and iterative algorithm is very effective here. 
Solution need not be unique depending in general on the 
starting point of iterations. Consequently, MaxLik es- 
timation provides a family of extremum states not dis- 
tinguished by the given measurement. Averaging over 
this family enhances the uncertainty of state prediction 
confirming the conjecture formulated in The analysis 
of realistic measurement supporting this interpretation is 
given in ||]. 

Relations (||) and (||) may show, how closely the given 
state approaches the extremum one. The relative entropy 
(normalized log likelihood) 



K{p/f) = -Y^fdnf>0 



(8) 



provides the difference between absolute minimum and 
estimated result. Its value may be expressed in % of the 
entropy S = — fi In fi . Similarly, the "experimentally 
achieved resolution of identity" R could always be com- 
pared with the identity operator. 

The MaxLik quantum state reconstruction possess 
very clear geometrical interpretation as normalization of 
incompatible observations. Indeed, the rays in Hilbert 
space are given up to the multiplicative factors. The 
renormalized projectors {y^l \y[) = \fjifpii\yi) fulfiU 



the relation analogous to (|l|) as 

kv\\p\v[) = n- 



(9) 



Moreover the operator R characterizes the overlapping 
of rays in analogy with tomography in medicine. The 
X-rays overlap in the scanned region illuminating larger 
part of the body. In the language of quantum theory 
the linear envelope of the detected projectors represents 
an analogy of the whole irradiated space U = {ly'i)}, a 
subspace of the full Hilbert space 7i. Denote formally 
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the orthogonal subspace of overlapping of projectors as 
O, \xk) being its orthogonal basis. Each projector can be 
then decomposed as jy-) = Yl^ki^klvDlxk) + \Zi), where 
\Zi) are orthogonal to subspace O, {xk\Zi) — 0. The sum 
of scalar products Z — J2i{Zi\Zi) characterizes the part 
of the projectors outside the orthogonal subspace O. The 
basis \xk) will be chosen in order to minimize Z under 
the condition of normalization. The optimum basis in O 
is spanned by the eigenstates diagonalizing the sum of 
projectors achieved by the realistic measurement 



\Xk) 



(10) 



Hence MaxLik reconstruction (Q) may be easily inter- 
preted in the language of quantum (Q) and geometri- 
cal (|o|) considerations. The operator R = J^ilvdiVil 
characterizes the overlapping of projectors. The recon- 
struction should be done in a subspace Orec spanned 
by the eigenstates \xk) with the degenerate eigenvalue 
Afc = 1. The density matrix is spanned in the subspace 
where R equals identity operator. The quantum postu- 
late (|l|) predicting the statistics of the outcome provided 
that quantum state is known should be modified to for- 
mally analogous relation (^ provided that an inversion 
problem has to be solved. However, the later problem is 
nonlinear since the normalization of projectors depends 
on the state itself. Any reconstruction beyond the sub- 
space Orec is rather a random guess since not enough 
information is available. In medicine, it would corre- 
spond to the "observation" of the head while stomach 
has been scanned, for example. The subspaces are re- 
lated as Orec C O C U C H, but may coincide in some 
special cases with low degrees of freedom, as for example 
in the case of spin systems. Since only the decomposi- 
tion of identity matters, the formulation is common for 
both the sharp and unsharp observations. The MaxLik 
estimations simplifies considerably, provided that mea- 
sured projectors H^ commute. Observations are compat- 
ible and renormalization is not necessary. The subspace 
O may be approximately characterized in the common 
basis 1^) by diagonal elements of operator R for pa — fi 
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The function R{^) plays the role of Optical Transfer 
Function in Fourier optics characterizing the fidelity of 
the information about variable ^ contained in the mea- 
sured data. 

The theory is illustrated on numerical simulations of 
the quantum tomography experiments of the type |^,|lO| . 
For quantum tomography, the projectors are given by 
rotated quadrature states \yij) = l/^/^^\xi,9j) corre- 
sponding to the center of the coordinate-phase bins. 
For random-phase homodyning the commuting projec- 
tors are Hi = l/(27r) J^^ d9\xi,9){xi,9\. Data depicted 
in the Fig. |l| correspond to the tomography of squeezed 



vacuum state-an eigenstate of the operator b = a cosh r-|- 
e^'^sinh ra^ for r = l,ip = Tr/2. Scanning has been done 
at 12 phase cuts using 600 records at each phase position. 
The scanned intervals (—7, 7) in each cut are divided into 
100 bins. An ideal detection {rj — 1) is assumed. Solution 
of the nonlinear eq. (H) represents the key point of the 
reconstruction. An iteration procedure has been already 
applied to random-phase homodyning |^,^. General it- 
eration procedure in tomography will be dealt with else- 
where. For the purpose of illustration the likelihood is 
maximized numerically using MATLAB for pure states. 
Estimated complex amplitudes in number-state ba- 
sis (hollow bars) are compared with the true amplitudes 
(full bars) in the left panels of the Fig. |[ The typical 
oscillating nature is obvious here. The hollow bars in 
right panels of the Fig. || show the diagonal elements of 
the operator R in 28 dimensional (upper panel) and in 
25 dimensional subspaces. The data are insufficient for 
MaxLik reconstruction in the former case since the iden- 
tity operator is not recovered, but are sufficient in the 
later case. For comparison, the full bars show the de- 
composition of identity for random-phase homodyning 
|ll|,Pt, where the reconstruction can be done in about 10 
dimensional subspace only. The relative entropy of the 
true state and generated data is about K{p/ f) = 0.95% 
of the entropy of measured data S — 5.89. The MaxLik 
fitting provides the value K{p/ f) — 0.81% of the entropy. 
For comparison, a random guess is characterized by the 
value of several tens of % of the entropy S. Finally, the 
Fi g. 1^ shows the pr oper normalization of the projectors 
^ f{xi, 9j)/ p{xi, 9j), by which the completeness relation 
is fulfilled. Not registered projectors are missing here 
and events registered with low relative frequencies are 
renormalized. 

MaxLik procedure provides an effective method for sta- 
tistical inverting of quantum postulate ([l]). It fits the 
data better than deterministic schemes, but may provide 
a family of indistinguishable states as the result. Max- 
Lik algorithm is nonlinear and standard error analysis 
cannot be applied. Particularly, reconstruction may be 
accomplished on the subspace where measured projectors 
provide the resolution of identity. These issues should be 
taken into account since the state reconstruction has be- 
came to play an important role in many sophisticated 
detection techniques. 
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FIG. 1. Histograms of homodyne detection for squeezed 
vacuum. 

FIG. 2. Left panels show the real and imaginary parts 
of the true (full) and estimated (hollow) amplitudes of pure 
state in number state basis. Right panels show the diagonal 
elements of decomposition of operator R for random-phase ho- 
modyne detection (full) and homodyne tomography (hollow) . 
For homodyne tomography, data are insufficient for recon- 
struction in 28 dimensional subspace (upper panel), but are 
sufficient to recover the identity in 25 dimensional subspace 
(lower panel). 

FIG. 3. Renormalization of projectors in the nonorthogo- 
nal rotated quadrature state basis corresponding to successful 
reconstruction. 
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HISTOGRAM OF HOMODYNE DETECTION 
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